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Abstract 

A simple statistical model in terms of light-front kinematic variables is used to 
explain the nuclear EMC effect in the range x S [0.2, 0.7], which was constructed by 
us previously to calculate the parton distribution functions (PDFs) of the nucleon. 
Here, we treat the temperature T as a parameter of the atomic number A, and get 
reasonable results in agreement with the experimental data. Our results show that 
the larger A, the lower T thus the bigger volume V , and these features are consistent 
with other models. Moreover, we give the predictions of the quark distribution 
ratios, i.e., q^{x)/q^ (x), q^[x)/q^[x), and s^{x)/s^{x), and also the gluon ratio 
g^{x)/g^{x) for iron as an example. The predictions are different from those by 
other models, thus experiments aiming at measuring the parton ratios of antiquarks, 
strange quarks, and gluons can provide a discrimination of different models. 
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1 Introduction 



Deep inelastic scattering (DIS) process is an efficient tool to detect the nucleon 
structure, and the high-energy charged lepton can also reveal the parton dis- 
tributions in the nuclear environment. The nucleons in a nucleus were initially 
thought to be highly insensitive to their surroundings, and the only nuclear 
effect in DIS was believed to be Fermi motion at large x. However, in 1982, 
the European Muon Collaboration at CERN discovered that nucleons inside 
a nucleus have a remarkably different momentum configuration as expected, 
i.e., the structure function for a bound nucleon differs from that for an iso- 
lated one significantly [TP|5] . which was named the nuclear EMC effect (for 
reviews, see, e.g., Refs. |I|5ll6H71l8]|9][T0] ) . The existence of the unexpected effect 
opens a window onto the behavior of quarks at nuclear scale, thus it offers the 
opportunity to use nuclear targets for exploring aspects of quark dynamics. 

In order to account for the EMC effect, there have been many efforts and in- 
sights implemented in various models, e.g., the cluster model |lllll2fl3fl4lll5fl6l 
the pion excess model p^l[T71[T8][T9ll20]|2T] . the x-rescaling model [22ll23]l2l] . 
the Q^-rescaling model [25|26f27j . the nucleon swelling model p3][23] . and the 
deconfinement model [T2|28] . In some sense, most of these available models 
provide a fairly good description, instead of an explanation, to the phenom- 
ena. The cluster model explains the EMC effect by introducing different ways 
of packing the constituents into clusters of 3n (n > 2) quarks, to mimic the 
field acting on partons caused by the nuclear mediate. This was suggested by 
Pirner et al. [TT], and ffist applied to the EMC effect by Jaffe [12]. Later it 
was discussed in detail by Carlson and Havens [13]. Furthermore, Barshay and 
Rein [T3] provided a physical picture for this scenario and predicted specific 
dependence. In the pion excess model. Smith [17] suggested that iron con- 
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tains 6~12 more pions which carry ~5% of the momentum, and Ericson and 
Thomas [IB] showed that the EMC enhancement in the region x < 0.3 can be 
reproduced with the excess pions if an attractive force was assumed. Close, 
Jaffe et al. |25|26|27] first discovered that the difference between the structure 
functions for a bound nucleon and that for a free one can be ascribed to the 
change of the effective Q^, which was named the Q^-rescahng model, and it 
was explained that these distinctions result of the difference in the 

scale of confinement of the nucleon. 

On the other hand, due to the complicated non-perturbative effect, it is still 
difficult to calculate the parton distribution functions (PDFs) of the free nu- 
cleon absolutely from the first principle theory of the quantum chromodynam- 
ics (QCD) at present. Various models according to the spirit of QCD have been 
brought forward, therein statistical ones, providing intuitive appeal and physi- 
cal simplicity, have made amazing success [29|30f31|l32|l33|l34f35f36f37|l38|l39f40f41f42f43f44f45f46f47f48|^ 
Actually, as can be speculated, with partons bound in the wee volume of the 
nucleon, not only the dynamic, but also the statistical properties, for example, 
the Pauli exclusion principle, should have important effect on PDFs. Angelini 
and Pazzi [29|30] . as pioneers, found that the nucleon valence quark distri- 
bution has a thermodynamical behavior for x > 0.1 with the temperature 
decreasing for different Q^. Cleymans and Thews [3T|32|33j explored a statis- 
tical way to generate compatible PDFs. Mac and Ugaz [M] incorporated first 
order QCD correction, and afterwards Bhalerao et al. [55]|36f37f38ll39] intro- 
duced finite-size correction; they both referred to the infinite-momentum frame 
(IMF). Devanathan et al. |^41f42] proposed a thermodynamical bag model 
which evolves as a function of x, and the structure functions they got have cor- 
rect asymptotic behavior; in addition, they parametrized on T and exhibited 
the scaling behavior. Bourrely, Soffer, and Buccella [I3]l44ll45f46|l47|48|49ll50|51f52] 
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developed a new form of statistical parametrization, and by incorporating 
QCD evolution they got indeed remarkable PDFs. Otherwise, Zhang et al. [5^115^ 
constructed a model using the principle of balance without any free parameter, 
and the Gottfried sum they got is in surprisingly agreement with experiments. 
Later, Singh and Upadhyay [SS] extended this model to have spin considered 
with modifications. Alberg and Henley [57] tracked the detailed balance model 
for a hadron composed of quark and gluon Fock states and obtained parton 
distributions for the proton as well as the pion. Recently, Trevisan et al. [58] 
presented a statistical model with a confining potential and took gluon split- 
ting into account further. 

The statistical idea is also applied to the nuclear EMC effect. Angelini and 
Pazzi [60] introduced thermodynamical analysis to the EMC effect, and by 
utilizing the ratio of valence quark distributions at different temperatures and 
confinement volumes, they fit data well. Afterwards, Li and Peng [61] discussed 
the EMC effect using the Fermi-Dirac distribution for fermions and Bose- 
Einstein distribution for gluons. Further, Rozynek and Wilk [62] combined 
nuclear Fermi motion with statistical effect [53|54|55j to account for the rise 
at large x. 

In Ref. [63], we performed a pure statistical model in terms of light-front 
kinematic variables, without any arbitrary parameter or extra corrected term 
put by hand, and obtained the analytic PDFs of free nucleons, i.e., the proton 
and the neutron. The results imply that the statistical effect is important to 
some aspects of the nucleon structure, and the simpleness and intuition of 
the model encourage us to apply it to the EMC effect. Here, we introduce 
the temperature T as a parameter versus the atomic number A, and the fit 
to the experimental data in the region x G [0.2, 0.7] is rather reasonable. 
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illustrated in Fig. [T], which indicating that the statistical effect may be an 
important source to the EMC phenomenon. But in the high-x region, where 
Fermi motion dominates to the effect, which is not considered in our model, 
our fit is below the data as expected. The temperature we get is about 1~2 
MeV lower in bound nucleons than in free ones, and jointly the volume is 
bigger about 5%~10%. Moreover, we also give the predictions of the quark 
distribution ratios, i.e., q^{x)/q^{x), q^{x)/q^{x), and s^{x) / and the 
gluon ratio g^{x)/g^{x) for iron as an example. 

2 Statistical approach 

We assume that the nucleon is a thermal system in equilibrium, made up 
of free partons (quarks, anti-quarks, and gluons), and the quarks and anti- 
quarks satisfy the Fermi-Dirac distribution while the gluons obey the Bose- 
Einstein distribution. Instead of boosting the distribution functions to the 
IMF |34|35ll36f37f38ll39j , we transform them in terms of light-front kinematic 
variables in the nucleon rest frame, and get x-dependent PDFs analytically [52] 




with the upper sign for fermions (quarks, anti-quarks), and nether sign for 
bosons (gluons); Qf is the degree of color-spin degeneracy, which is 6 for quark 
(anti-quark) and 16 for gluon; /i/ is the corresponding chemical potential. 



5 



while for anti-quark /z^ = — /iq, and for gluon /ig = 0; x is the hght-front 
momentum fraction of the nucleon carried by the specific parton; M is the 
mass of the nucleon, and the value is taken as 938.27 MeV; and Ijhiz) is the 
polylogarithm function, defined as Li2(z) = HkLi /k"^- Note that the analytic 
expression above is different from those attained in the previous statistical 
models [n]l32]l3Mll35]l3M71l38]^^ 



The expression of the PDFs we derived is available for all the partons with 
mass mj, however, for simplicity, in Ref. [53], we mainly focused on the u, 
d flavors, which can be viewed massless as a quite good approximation. In 
this paper, we will take a further consideration of the s flavor, whose mass is 
around 100 MeV, so the constraint x > m'l/M'^ should be fulfilled [65]. For 
convenience of a universal manipulation to all partons, we have multiplied 
explicitly to the expression of the PDFs, z.e., Eq. ([1]), a factor of the step- 
function, 9{x — im?r/M'^). 



In practice, the PDFs in a certain system should be constrained with some 
conversation laws. For example, in the proton, they are 

uv = j [u{x) - u{x)] dx = 2 , (2) 

dy = J [d{x) — d{x)] dx = 1 , (3) 
fxf{x)dx = l, (4) 

where /, in Eq. (j4]), denotes to u, u, d, d, g if just the u, d flavors and the 
gluon are considered, and to u, u, d, d, s, s, g when the s flavor is appended. 
However, the contribution of the heavier flavors, the c quark for instance, is 
negligible [63] . 
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There are four unknown parameters T, V, jJ'd and three constraints, thus 
for a fixed T, the rest parameters V, fiu, /id can be determined uniquely, 
whereafter we can obtain all the PDFs. Worthy to mention that, including 
the s flavor will not introduce any extra parameter, due to /xg = /i^ = 0. 

After the PDFs of the proton have been determined, we can get the PDFs 
of the neutron as well, using the p-n isospin symmetry, i.e., = (P{x), 

= nP{x), vT'^x) = (F{x), cr(x) = vP{x), g'^{x) = gP{x), and s^{x) = 

= sP{x) = sP{x). 

Eqs. (l2])-(j4]) should be satisfied for free nucleons, as well as the ones immersing 
in the nuclear environment. Here we mainly assume that a nucleon under a 
different nuclear circumstance is equivalent to at a different temperature, and 
subsequently along with different V, fi^, and fid- 

Also note that, the PDFs and the nucleon structure function 

F2{x) = xY.e}f{x), (5) 
/ 

where e/ is the charge of the parton of flavor /, as well as the Gottfried sum 

o _ r' F'iix) - F-{x) ^ 1 2rir . 

Sg = / dx = - - - d{x} - u{x) dx , (6) 

Jo X 6 6 Jo ^ J 

are actually all T-dependent. 

Afterward, we perform the numerical calculation to determine the tempera- 
tures of the nucleons in different nuclear environment, first in deuteron, then 
in other nuclei. Since deuteron is a weakly bound system, the structure func- 
tion of it can be taken as a free isoscalar nucleon approximately, i.e., F^{x) = 
[Ff (x) + -F^(x)] /2. In practice, we figure out T^, i.e., the temperature of nu- 
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cleon in deuteron by using Sg(T^) = Sq^, and then fit F^{T^, x) / (T^ , x) 
to available experimental data in the region x G [0.2, 0.7] where the EMC ef- 
fect dominates, thus the parameters of the nucleon in the nucleus with atomic 
number A can be determined. 

3 Results and Discussions 

Numerical calculation is performed under three different conditions below 

• just considering the u, d flavors and the gluon, 

• including the s flavor and taking the mass of s quark = 130 MeV, 

• including the s flavor and taking = 70 MeV; 

here 130 MeV and 70 MeV are respectively the upper and lower limit of the 
mass of s quark suggested in PDG2008 [66] . 

In practice, as mentioned above, we figure out the temperature of the nucleon 
in deuterium firstly, at which the Gottfried sum Sq equals to the exper- 
imental value 5*^^^ = 0.235 ± 0.026 [67f68] . According to the three different 
conditions listed above, the results are a bit different 

• = 46.84 MeV when only u, d and g are considered, 

• = 42.53 MeV when u, d, s and g are considered with = 130 MeV, 

• = 40.08 MeV when u, d, s and g are considered with = 70 MeV. 

Then we fit the analytic expression F.f^(T^,x)/F^(T^ ,x) to various experi- 
mental data from Ref. |6l] in the region x G [0.2, 0.7] under the constraints 
aforementioned, and we attain the corresponding temperatures of the nucleons 
inside various nuclei. The fitting curves of some nuclei are shown in Fig. [H 
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We can see that, when x decreases from 0.2 to about 0.1, our curves predict 
the right trend and fit to the experimental data roughly; but confiicted with 
the data, as x decreases further from 0.1, our curves still go up and stand as 
a constant slightly above unity when x approaches 0. In the high- a; direction, 
the experimental points reveal an apparent rise from about x = 0.7~ 0.8 as x 
increases, however, our curves depart from them and monotonously go down, 
caused by the neglect of Fermi motion of the nucleons. 

As generally recognized now, the fact that =^ F.f is referred to with differ- 
ent terms: the behavior of F.^{x)/F^{x) in the region 0.2 < x < 0.7 is often 
referred to as the "(special) EMC effect", where the ratio goes down straight 
to lower than unity; in the region x < 0.05^0.1, the ratio is smaller than 
unity, named "shadowing", and in the range 0.1 < x < 0.2, the ratio shows an 
enhancement and becomes a little larger than unity, called "anti-shadowing" ; 
further, in the region x > 0.7^0.8, the ratio has a precipitous rise and grows to 
above unity, where is the mentioned "Fermi motion" region. Here, our model 
describes the "EMC effect" reasonably but fails in other regions, like most of 
the other models not considering the effect of "Fermi motion" and "shadowing 
effect". 

Once the temperatures of the nuclei are determined, the other parameters are 
also fixed simultaneously concerning Eqs. and the values are listed in 

Tables [H [21 [3] respectively according to the three different conditions noted 
above. For an intuitive view, we present the physical quantities T, So, V, r, 
Hu, and Hd in Fig. [21 where r = (SV/AttY^'^ represents the radius of nucleon 
from the model. We can see clearly that the nucleus with larger A has a 
lower temperature T, jointly a larger volume V and smaller quark chemical 
potentials fid- 
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The features derived above have an self-consistent explanation in a physical 
picture. When A increases, the volume, in which the partons move, becomes 
larger, thus the average kinetic energy of the partons decreases, due to the 
uncertainty principle. This agrees with the corresponding lower T property. On 
the other hand, as V increases, the number densities of quarks decrease, thus 
the chemical potentials, which are monotonic functions of the corresponding 
number densities, become lower. 

Our result that the volume of the nucleon in nuclei becomes larger is quali- 
tatively consistent with other models, such as the Q^-rescaling model and the 
nucleon swelling model, though the starting point of our model is different. 
The result we derive, - = 1~2 MeV, is very close to 3 ± 1 MeV given 
in Ref. [60], and our enhancement in V'^/V^ is about half of that presented 
in Ref. [61]. Worthy to note that, including the s flavor and taking differ- 
ent masses of it lead to some slightly difference in results, so the s flavor is 
considered as a modification here. 

In addition, as parameters are determined, we can give explicitly the predic- 
tions of PDFs of the nucleons inside different nuclei. The ratios of the PDFs 
of iron to deuterium are depicted in Fig. [3] for example. Note that, the distri- 
bution of d flavor, not presented in Fig. [31 is just the same as that of u flavor, 
for that all nuclei here, including iron, have been adjusted to be isospin scalar. 
The solid normal, thick, and thin curves correspond respectively to our results 
in the three conditions mentioned above, and they almost overlap in each sub 
figure, though the parameters of them are a bit different (see Tables [1], [21 [3]). 
The sub figure of s"^/ in Fig. [3] lacks the normal solid curve, for it is the 
condition where no s flavor is included. 
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The ratios in Fig. [3] are only predicted in the range x G [0.2, 0.7] since we 
only fit to the data in this special region and do not include other effect. We 
can see that the ratios are monotonous of x, indicating more low-x momen- 
tum configuration in iron than in deuterium, which is in consistency with our 
intuition of a lower temperature in iron. 

And also, we present the fractions of momentum x and structure function F2 
carried by various partons under the three conditions respectively in Table HI 
From the results, we can see that as getting smaller, s quarks take up 
more fractions of momentum and structure function, and the contributions 
from which are somehow considerable. The case with no s quarks can be 
considered as a limit for 00. 

Together with our results for iron, the predicted PDFs of other models, i.e., 
the cluster model, the pion excess model, and the Q^-rescaling model, are also 
illustrated in Fig. O 

For the cluster model, here we consider 6-quark cluster only, and employ the 
PDFs given in Ref . [16] . Also we assume that the probability to form a 6-quark 
cluster in iron equals to 0.3 [13]. 

Concerning the pion model, we adopt the CTEQ6M parametrization [69] for 
nucleon and the RMS parametrization [2U] for pion, then the PDFs in nuclei 
read [T7IIT81I70] 

<if{^) = f ^f^y)^n-) + f ^/^(y)gf (-) , (7) 

Jx y y Jx y y 

where qf{x) and qf{x) are respectively the parton distributions for the free 
pion and nucleon, and we make use of the "toy model" in Ref. [21] to give the 
probability fn{y) and f^{y) of finding pionic and normal nucleonic content 
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respectively in iron. 



And as to the Q^-rescaling model, we again take the CTEQ6M parametriza- 
tion [Sn] and take C = 1-83 [27] in the PDFs below 



g^(x,g^) = g^(x,e(g^)g^). 



As can be seen in Fig. [3], the predicted PDFs differ significantly from model to 
model. And to distinguish various models and look into the immanent reason 
of the EMC effect, we suggest more experiments to identify the PDFs in nuclei, 
especially for anti-quarks, the strange quark, and the gluon. The dimuon yield 
in Drell-Yan process [7T|72] can detect the sea content in nuclei. Ref. [70] also 
advised that the semi-inclusive hadron productions in DIS are sensitive to the 
sea quark content. Further, the process of charmed quarks production in DIS 
can probe the gluon constituent via the photon-gluon fusion mechanism, and 
the strange quark content can be detected in A-K process. 



4 Summary 



In this paper, we utilize a concise statistical model to mimic the nuclear EMC 
effect. The statistical model was used previously to calculate the parton dis- 
tribution functions (PDFs) of the nucleon in terms of light-front kinematic 
variables with no arbitrary parameter or extra corrected term [63]. Here, we 
treat the nucleon temperature T as a parameter of the atomic number A and 
find that the nuclear effect can be explained as a shift of T — as A increases, 
T gets lower and consequently the volume V becomes larger; the larger A, the 
more significant influence. These features are consistent with other models. 
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for example, the Q^-rescaling model and the nucleon swelling model. 

Further, we present the predictions of the quark distribution ratios, i.e., 
q^{x)/q^{x), and s"^ (x) / (x) , and also the gluon ratio g^{x)/g^{x) for iron 
as an example. These predictions are rather different from those of other avail- 
able models. Further experiments are expected to provide more information 
of the PDFs in nuclei, especially for anti-quarks, the strange quark, and the 
gluon, then we can test various models better. 

Though there is only one free parameter, i.e., the nucleon temperature T, in 
our model, the descriptions to the experimental data are reasonable in the 
EMC domain x G [0.2, 0.7]. We suggest that the statistical property may be 
a considerable source to the nuclear EMC effect, and expect to have more 
relevant experiments to test the predictions. 

Acknowledgment 

This work is partially supported by National Natural Science Foundation of 
China (No. 10721063) and by the Key Grant Project of Chinese Ministry of 
Education (No. 305001). It is also supported by Hui-Chun Chin and Tsung- 
Dao Lee Chinese Undergraduate Research Endowment (Chun-Tsung Endow- 
ment) at Peking University, and by National Fund for Fostering Talents of 
Basic Science (Nos. J0630311, J0730316). 

References 

[1] J.J. Aubert et al, Phys. Lett. B 123 (1983) 275. 



13 



[2] J.J. Aubert et al, Nucl. Phys. B 293 (1987) 740. 

[3] A. Bodek et al, Phys. Rev. Lett. 50 (1983) 1431. 

[4] E.L. Berger and F. Coester, Ann. Rev. Nucl. Part. Sci. 37 (1987) 463. 

[5] L.L. Frankfurt and M.I. Strikman, Phys. Rep. 160 (1988) 235. 

[6] V. Barone and E. Predazzi, Ann. Phys. -Paris 12 (1987) 525. 

[7] M. Arneodo, Phys. Rep. 240 (1994) 301. 

[8] D.F. Geesaman, K. Saito, and A.W. Thomas, Ann. Rev. Nucl. Part. Sci. 45 
(1995) 337. 

[9] G. Filler and W. Weise, Phys. Rep. 330 (2000) 1. 

[10] P.R. Norton, Rep. Prog. Phys. 66 (2003) 1253. 

[11] H.J. Pirner and J.P. Vary, Phys. Rev. Lett. 46 (1981) 1376. 

[12] R.L. Jaffe, Phys. Rev. Lett. 50 (1983) 228. 

[13] C.E. Carlson and T.J. Havens, Phys. Rev. Lett. 51 (1983) 261. 

[14] S. Barshay and D. Rein, Z. Phys. C 46 (1990) 215. 

[15] J. Bias de Deus, M. Pimenta, and J. Varela, Z. Phys. C 26 (1984) 109. 

[16] U. Sukhatme, G. Wilk, and K.E. Lassila, Z. Phys. C 53 (1992) 439. 

[17] C.H.L. Smith, Phys. Lett. B 128 (1983) 107. 

[18] M. Ericson and A.W. Thomas, Phys. Lett. B 128 (1983) 112. 

[19] B.L. Friman, V.R. Pandharipande, and R.B. Wiringa, Phys. Rev. Lett. 51 
(1983) 763. 

[20] P.J. Sutton, A.D. Martin, R.G. Roberts, and W.J. Stirling, Phys. Rev. D 45 
(1992) 2349. 

14 



[21] E.L. Berger, F. Coester, and R.B. Wiringa, Phys. Rev. D 29 (1984) 398. 

[22] M. Staszel, J. Rozynek, and G. Wilk, Phys. Rev. D 29 (1984) 2638. 

[23] A.W. Hendry, D.B. Lichtenberg, and E. Predazzi, Nuovo Cim. A 92 (1986) 
427. 

[24] S.V. Akulinichev, S. Shlomo, S.A. Kulagin, and G.M. Vagradov, Phys. Rev. 
Lett. 55 (1985) 2239. 

[25] F.E. Close, R.G. Roberts, and G.G. Ross, Phys. Lett. B 129 (1983) 346. 

[26] R.L. Jaffe, F.E. Close, R.G. Roberts, and G.G. Ross, Phys. Lett. B 134 (1984) 
449. 

[27] F.E. Close, R.L. Jaffe, R.G. Roberts, and G.G. Ross, Phys. Rev. D 31 (1985) 
1004. 

[28] T. Kawabe, Lett. Nuovo Cim. 44 (1985) 159. 

[29] C. Angelini and R. Pazzi, Phys. Lett. B 113 (1982) 343. 

[30] C. Angelini and R. Pazzi, Phys. Lett. B 135 (1984) 473. 

[31] J. Cleymans, Presented at International Workshop on Quarks, Gluons and 
Hadronic Matter, Cape Town, South Africa, 16-20 February 1987, published 
in Cape Town Workshop 1987, 0447 (QCD161:C375:1987). 

[32] J. Cleymans and R.L. Thews, Z. Phys. C 37 (1988) 315. 



[33] J. Cleymans, L Dadic, and J. Joubert, arXiv:hep-ph/9307352 



[34] E. Mac and E. Ugaz, Z. Phys. C 43 (1989) 655. 



[35] R.S. Bhalerao, arXiv:hep-ph/9506367 



[36] R.S. Bhalerao, Phys. Lett. B 380 (1996) 1. 
[37] R.S. Bhalerao, Nucl. Phys. A 680 (2000) 62. 

15 



[38] R.S. Bhalerao, Phys. Rev. C 63 (2001) 025208. 

[39] R.S. Bhalerao, N.G. Kelkar, and B. Ram, Phys. Lett. B 476 (2000) 285. 

[40] K. Ganesamurthy, V. Devanathan, and M. Rajasekaran, Z. Phys. C 52 (1991) 
589. 

[41] V. Devanathan, S. Karthiyayini, and K. Ganesamurthy, Mod. Phys. Lett. A 9 
(1994) 3455. 

[42] V. Devanathan and J.S. McCarthy, Mod. Phys. Lett. A 11 (1996) 147. 

[43] C. Bourrely et al, Z. Phys. C 62 (1994) 431. 

[44] C. Bourrely and J. Soffer, Phys. Rev. D 51 (1995) 2108. 

[45] C. Bourrely and J. Soffer, Nucl. Phys. B 445 (1995) 341. 

[46] C. Bourrely and J. Soffer, Phys. Rev. D 68 (2003) 014003. 

[47] C. Bourrely, J. Soffer, and F. Buccella, Eur. Phys. J. C 23 (2002) 487. 

[48] C. Bourrely, J. Soffer, and F. Buccella, Eur. Phys. J. C 41 (2005) 327. 

[49] C. Bourrely, J. Soffer, and F. Buccella, Mod. Phys. Lett. A 21 (2006) 143. 

[50] C. Bourrely, J. Soffer, and F. Bucceha, Phys. Lett. B 648 (2007) 39. 

[51] C. Bourrely, J. Soffer, and F. Buccella, Mod. Phys. Lett. A 18 (2003) 771. 

[52] J. Soffer, Nucl. Phys. A 755 (2005) 361. 

[53] Y.-J. Zhang, B. Zhang, and B.-Q. Ma, Phys. Lett. B 523 (2001) 260. 
[54] Y.-J. Zhang, B.-S. Zou, and L.-M. Yang, Phys. Lett. B 528 (2002) 228. 
[55] Y.-J. Zhang, W.-Z. Deng, and B.-Q. Ma, Phys. Rev. D 65 (2002) 114005. 
[56] J.P. Singh and A. Upadhyay, J. Phys. G 30 (2004) 881. 
[57] M. Alberg and E.M. Henley, Phys. Lett. B 611 (2005) 111. 

16 



[58] L.A. Trevisan, C. Mirez, T. Frederico, and L. Tomio, Eur. Phys. J. C 56 (2008) 
221. 

[59] R.P. Bickerstaff and J.T. Londergan, Phys. Rev. D 42 (1990) 3621. 
[60] C. Angelini and R. Pazzi, Phys. Lett. B 154 (1985) 328. 

[61] Z. Li and H. Peng, High Energy Phys. Nucl. Phys. (in Chinese) 12 (1988) 55. 

[62] J. Rozynek and G. Wilk, Nucl. Phys. A 721 (2003) 388. 

[63] Y. Zhang, L. Shao, and B.-Q. Ma, Phys. Lett. B 671 (2009) 30. 

[64] J. Gomez et al, Phys. Rev. D 49 (1994) 4348. 

[65] P. Zavada, Phys. Rev. D 55 (1997) 4290. 

[66] Particle Data Group, Phys. Lett. B 667 (2008) 1. 

[67] New Muon Collaboration, P. Amaudruz et al, Phys. Rev. Lett. 66 (1991) 2712. 

[68] New Muon Collaboration, M. Arneodo et al, Phys. Rev. D 50 (1994) Rl. 

[69] J. Pumplin et al, JHEP 0207 (2002) 012. 

[70] B. Lu and B.-Q. Ma, Phys. Rev. C 74 (2006) 055202. 

[71] D.M. Aide et al, Phys. Rev. Lett. 64 (1990) 2479. 

[72] G. T. Garvey and J. C. Peng, Prog. Part. Nucl. Phys. 47 (2001) 203 
jarXiv:nucl-ex/0109010] . 



17 



List of Figures 

• Fig.l: The ratios of the structure function F^jF^- The sohd, dashed, 
and dotted curves correspond to three conditions: a) just considering the 
M, d flavors and the gluon, b) considering the u, ci, s flavors and the gluon 
with rtis = 70 MeV, c) considering the u, d, s flavors and the gluon with 
rris = 130 MeV, respectively. The experimental data are from Ref. [M] . 

• Fig. 2: The parameters of the nucleons in nuclei — T, Sq, V^/V^, r^/r^, 
fiu, and fid versus the atomic number A. T, fiu, and /i^ are in units of 
MeV. The squared, circular, and triangular points correspond to the three 
conditions: a) just considering the u, d flavors and the gluon, b) considering 
the u, d, s flavors and the gluon with rris = 70 MeV, c) considering the u, 
d, s flavors and the gluon with = 130 MeV, respectively. 

• Fig. 3: The ratios of the PDFs of iron to deuterium. The solid normal, 
thick, and thin curves are our results, corresponding respectively to the three 
conditions: a) just considering the u, d flavors and the gluon, b) considering 
the u, d, s flavors and the gluon with rris = 70 MeV, c) considering the u, 
d, s flavors and the gluon with rris = 130 MeV. The dashed, dotted, and 
dash-dotted curves correspond to the predictions of the cluster model, the 
Q^-rescaling model, and the pion excess model, respectively. 



18 



List of Tables 



• Table 1: Considering the m, d flavors and the gluon. 

• Table 2: Considering the u, d, s flavors and the gluon (m^ = 130 MeV). 

• Table 3: Considering the u, d, s flavors and the gluon {rris = 70 MeV). 

• Table 4: The fractions of x and F2 carried by u, u, d, d, g, s{s) in 
and ^^Fe respectively under the three conditions: a) just considering the u, 
d flavors and the gluon, b) considering the u, d, s flavors and the gluon 
with rris = 70 MeV, c) considering the u, d, s flavors and the gluon with 

= 130 MeV. 



19 




Fig. 1. The ratios of the structure function jF^ ■ The sohd, dashed, and dotted 
curves correspond to three conditions: a) just considering the tt, d flavors and the 
gluon, b) considering the n, s flavors and the gluon with = 70 MeV, c) 
considering the u, d, s flavors and the gluon with = 130 MeV, respectively. The 
experimental data are from Ref. |64j . 
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Fig. 2. The parameters of the nucleons in nuclei — T, Sq, V'^/V^ , j^D ^ ^^^^ 
Hd versus the atomic number A. T, and Hd are in units of MeV. The squared, 
circular, and triangular points correspond to the three conditions: a) just considering 
the M, d flavors and the gluon, b) considering the ti, d, s flavors and the gluon with 
rus = 70 MeV, c) considering the u, d, s flavors and the gluon with = 130 MeV, 
respectively. 
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Fig. 3. The ratios of the PDFs of iron to deuterium. The solid normal, thick, and 
thin curves are our results, corresponding respectively to the three conditions: a) 
just considering the u, d flavors and the gluon, b) considering the u, d, s flavors 
and the gluon with nig = 70 MeV, c) considering the u, d, s flavors and the gluon 
with nis = 130 MeV. The dashed, dotted, and dash-dotted curves correspond to 
the predictions of the cluster model, the Q^-rescaling model, and the pion excess 
model, respectively. 
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Table 1 



Considering the u, d flavors and the gluon. 





2D 9Be ^^C ^7^1 40Ca sepg lOT^g wtj^^ 


T (MeV) 


46.84 46.00 45.97 45.76 45.64 45.64 45.25 45.29 45.01 


/in (MeV) 


63.21 60.67 60.58 59.96 59.60 59.60 58.46 58.58 57.77 


/id (MeV) 


35.39 33.83 33.77 33.39 33.17 33.17 32.48 32.55 32.06 




1.000 1.088 1.091 1.114 1.128 1.128 1.174 1.169 1.203 




1.000 1.028 1.029 1.037 1.041 1.041 1.055 1.053 1.063 


Sg 


0.235 0.232 0.231 0.231 0.230 0.230 0.228 0.229 0.227 
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Table 2 



Considering the u, d, s flavors and the gluon (m^ = 130 MeV). 





2D ^He 12^ 27^1 40Ca sep^ lOT^g 197^^ 


T (MeV) 


42.53 41.87 41.84 41.69 41.58 41.59 41.27 41.30 41.07 


/in (MeV) 


57.39 55.31 55.21 54.75 54.41 54.44 53.46 53.55 52.86 


/id (MeV) 


32.13 30.85 30.79 30.50 30.29 30.31 29.71 29.77 29.34 




1.000 1.078 1.081 1.100 1.114 1.112 1.154 1.150 1.181 




1.000 1.025 1.026 1.032 1.037 1.036 1.049 1.048 1.057 


Sg 


0.235 0.232 0.232 0.231 0.230 0.230 0.229 0.229 0.228 
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Table 3 



Considering the u, d, s flavors and the gluon (m^ = 70 MeV). 





2D ^He 12(3 40ca sep^ lOT^g 197^^ 


T (MeV) 


40.08 39.52 39.50 39.37 39.28 39.28 39.01 39.04 38.84 


/in (MeV) 


54.08 52.39 52.33 51.95 51.68 51.68 50.88 50.97 50.39 


/id (MeV) 


30.28 29.24 29.20 28.97 28.80 28.80 28.32 28.37 28.01 




1.000 1.067 1.070 1.086 1.098 1.098 1.133 1.129 1.157 




1.000 1.022 1.023 1.028 1.032 1.032 1.043 1.041 1.050 


Sg 


0.235 0.232 0.232 0.232 0.231 0.231 0.230 0.230 0.229 
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Table 4 

The fractions of x and F2 carried by u, u, d, d, g, s{s) in and ^^Fe respectively 
under the three conditions: a) just considering the u, d flavors and the gluon, b) 
considering the u, d, s flavors and the gluon with nis = 70 MeV, c) considering the 



u, d, s flavors and the gluon with lUg = 130 MeV. 



The fraction 


u u d d g s{s) 


2D 

X 

no s quark 

F2 


0.295 0.041 0.295 0.041 0.328 
0.289 0.043 0.289 0.043 0.336 
703 097 1 76 024 

W.IWO W.WC/I W.±lvJ W.WiiiT: 

0.696 0.104 0.174 0.026 


2D 

X 

56pg 

m, = 70 MeV 

2D 

56Fe 


0.253 0.035 0.253 0.035 0.281 0.072 
0.249 0.037 0.249 0.037 0.286 0.072 
0.640 0.088 0.160 0.022 0.045 
0.634 0.093 0.158 0.023 0.046 


2D 

X 

= 130 MeV 

2D 
56pg 


0.268 0.037 0.268 0.037 0.298 0.046 
0.263 0.039 0.263 0.039 0.305 0.045 
0.664 0.091 0.166 0.023 0.028 
0.657 0.098 0.164 0.025 0.028 
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